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Abstract 

Decompositions of the unitary group U(n) are useful tools in quantum information the- 



■ ory as they allow one to decompose unitary evolutions into local evolutions and evolutions 

causing entanglement. Several recursive decompositions have been proposed in the literature 
£f} , to express unitary operators as products of simple operators with properties relevant in en- 

^\ • tanglement dynamics. In this paper, using the concept of grading of a Lie algebra, we cast 

these decompositions in a unifying scheme and show how new recursive decompositions can 
be obtained. In particular, we propose a new recursive decomposition of the unitary operator 
on TV qubits, and we give a numerical example. 
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5 '. 1 Introduction 

q-( Decompositions of the unitary Lie group U(n) serve to factorize any element Xf 6 U(n) as a 

product Xf — XiX 2 ■ ■ ■ X m , where X\, . . . , X m are (elementary) factors in U (n). There are several 
reasons to study such decompositions for unitary evolutions in quantum mechanics. They allow one 
to analyze the dynamics of a quantum system in terms of simpler, possibly meaningful, factors. In 
particular, for multipartite systems they allow the identification of the local and entangling parts of 
a given evolution. In this context one can study entanglement dynamics pQ , [5] , [H] ■ From a more 
practical point of view, they allow one to decompose the task of designing a given evolution, such as 
a quantum gate, into simpler, readily available dynamics (cf., e.g., [H]). In particular, in multipartite 
systems few entangling evolutions are typically available. Lie group decompositions are also useful 
in control problems [3] and in the solution of some algebraic problems of interest in quantum 
information [S]. For these reasons several decompositions have been introduced in recent years 
[I]) 0) [3]> E], U\- I n DP' 12' a decomposition called the Concurrence Canonical Decomposition 
(CCD) was studied in the context of entanglement theory. The CCD is a way to decompose 
every unitary evolution on N q-bits into a part that does not modify the concurrence on the iV 
q-bits, and a part that does. It is a Cartan decomposition in that it corresponds to a symmetric 
space of SU(2 N ) [B]. In [3] the CCD was further studied and generalized to multipartite systems of 
arbitrary dimensions. The resulting decomposition was called an Odd-Even Decomposition ( OED). 
The OED is a decomposition of unitary evolutions on multipartite systems constructed in terms 
of decompositions on the single subsystems. Recursive decompositions such as the ones in [5] and 
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[7J recursively apply the Cartan decomposition theorem in order to decompose the factors into 
simpler ones. 

The present paper is devoted to recursive decompositions. Using the relation between Cartan 
decompositions of Lie algebras and Lie algebra gradings, we show that the recursive decompositions 
of [5] and [Jj are a special case of a general scheme from which several other recursive decompositions 
can be obtained. 

The paper is organized as follows. Most of the content of section[2]is background material con- 
cerning the basic concepts of Cartan decompositions of Lie groups and algebras, with particular 
emphasis on decompositions of U(n). We also describe the main ingredients of the CCD decomposi- 
tion of [1] and [2]; the OED decomposition of [4]; and the recursive decompositions of Khaneja and 
Glaser [7J, and D'Alessandro and Romano [5J. One extension of the procedure used for the OED 
decomposition is presented in Theorem 1 2. II In section[31 we describe gradings of Lie algebras, and 
establish a link between gradings and recursive decompositions. This gives a general method to 
develop recursive decompositions of U(n). We show in section [¥] how the recursive decompositions 
of [5] and [7] are special cases of this general procedure, and how new recursive decompositions can 
be obtained. In section O we give a numerical example illustrating the calculation of the recursive 
decompositions described in section |4] Some concluding remarks are presented in section [6l 



2 Cartan decompositions of the unitary group 

2.1 Cartan decompositions of a Lie Algebra 

A Cartan decomposition of a semisimple Lie algebra £ is a vector space decomposition 

c = k®v, (i) 

where the subspaces K, and V satisfy the commutation relations 

[K,K\QK, [K,,V]QV, [P,V]CK,. 

The pair (JC,P) is called a Cartan pair of C. In particular, K, is closed under the Lie bracket, 
and is therefore a Lie subalgebra of C. A Cartan decomposition of a Lie algebra C induces a 
decomposition of the connected Lie group associated to £, which we denote by e . In particular, 
every element L of e c can be written as 

L = KP, (2) 

where K € e K and P is the exponential of an element in V . Since [P \V\ C IC, any Lie subalgebra 
contained in V is necessarily Abelian. A maximal Abelian subalgebra TL contained in V is called a 
Cartan subalgebra, and the common dimension of all the maximal Abelian subalgebras TL is called 
the rank of the decomposition. Indeed, although the Cartan subalgebra is not unique, it may be 
shown that two Cartan subalgebras Tt and TLi are conjugate via an element of e . This means 
that there exists S S e such that TL = Adg(TLi). Here Ads denotes the adjoint map defined as 
Ad s (H) := SHS^ for H e C. 

Let TL be a Cartan subalgebra of C. One can prove that 

V= (J Ad s (H), 

and therefore 

exp(P) = |J Ad s (e n ). 

It follows that P in (J2} has the form P = SAS^, with S £ e K and A E e 
L £ e c can be written as 

L = K t AK 2 , 



Hence, from ([2]), each 
(3) 



3 



where K\, K2 € e K and A e e n . This decomposition is known as the KAK decomposition of the 
Lie group e c . 

Cartan classified all the Cartan decompositions of the classical Lie algebras [BJ. In particular, 
up to conjugacy, there exist three types of Cartan decomposition of the special unitary Lie algebra 
su(n) , the Lie algebra of skew-Hermitian matrices with zero trace. The decompositions are classified 
as AI, All, and AIII. 

A decomposition of type AI is the Cartan decomposition of su(n) into purely real and purely 
imaginary matrices, i.e., 

Su(n) = so(n) © so(n)- 1 -. (4) 

The orthogonality is given by the inner product (A , B) = Tr(ABt) where A,Be su(n). The 
diagonal matrices in so(n) span a maximal Abelian subalgebra, so the rank of the decomposition 
is rt — 1. 

A decomposition of type All is of the form 

su(2n) = sp(n) © sp(n)- L , (5) 

where sp(n) is the Lie algebra of symplectic matrices, namely the subalgebra of su(2n) of matrices 
A satisfying 

AJ + JA T = 0, 



in which J is the 2n x In matrix 



J 



l n 

-In 



Here and in the rest of this paper, we denote by l n the n x n identity matrix. The rank of the 
decomposition All is again n — 1. 

A decomposition of type AIII is defined in terms of two positive integers p and q with p + q = n. 
The decomposition is 

su{n) := JC®V, (6) 
where JC is spanned by block diagonal matrices 

with X pxp and Y qxq skew-Hermitian and Tr(Ap X p) + Tr(y gxg ) = 0. The rank of this decomposition 
is min{p, q}. 

Each Cartan decomposition of su(n) is conjugate to one of the decompositions of type AI, All, 
and AIII. In other words, if su(n) = JC © V is a Cartan decomposition of su(n), there exists a 
unitary matrix T such that fiu(n) = © TVT^ is in one of the forms AI, All and AIII. These 

decompositions can be expressed in forms of interest in various contexts, for example with matrices 
expressed as tensor products of operators on single subsystems in a multipartite quantum system. 

In the following, we shall find it convenient to extend these decompositions to decompositions 
of u(n) = su(n) © span{zl„}, the Lie algebra of U(n). Consider a Cartan decomposition of the 
special unitary Lie algebra su(n) of type either AI or AIL Since the identity matrix 1„ commutes 
with each element of su(n) , the Cartan decompositions of su(n) of types AI ((4| and All (|5|) can be 
naturally extended to decompositions of u(n) by replacing V with ^©spanli 1„}. We also denote 
these decompositions of types AI and AIL In both Cartan decompositions, the rank becomes n. 
For decompositions of type AIII, we find it convenient to include span{i 1„} in the Lie algebra 
part, and replace JC with JC © span{i l n }, so as to lift the restriction Tr(A pxp ) + Tr(Y qxq ) = in 
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2.2 Cartan decompositions for multipartite quantum systems; CCD and 
OED 

For a multipartite quantum system with N subsystems of dimensions ni,...,njv, the set of pos- 
sible Hamiltonians is the Jordan algebra iu(niri2 • ■ ■ njv) of ri\ni • ■ ■ njsr x nxn% ■ ■ ■ Hermitian 
matrices. The Lie algebra associated to the dynamics is u{n\ri2 • ■ ■ Ujv). Cartan decompositions of 
u(n\n2 ■ ■ ■ njsr) result in decompositions of the corresponding unitary group of quantum evolutions 
U(nin 2 ■ ■ ■ n N ). 

The Concurrence Canonical Decomposition (CCD) was studied in pQ [2] as a means of decom- 
posing the dynamics of iV two level systems, into one factor which preserves the concurrence of 
the density matrix, and one factor which does not. It is constructed as follows: 

Recall that the Pauli matrices 

^ = (i o)' ^ = (° o*)' °"* = (o -l)' 

together with the 2x2 identity matrix I2, form a basis of the Jordan algebra iu(2). An orthogonal 
basis of u(2 N ) is given by the tensor products of the form i<j\ ® ■ ■ ■ ® cat, where dj = <Jx.y,z or 
dj = I2 for all 1 < j < N. Let us denote by iX a and iX e the respective subspaces of u(2 JV ) spanned 
by elements of the form i<j\ (E> ■ ■ • <E> cjv with an odd or even number of factors <jj given by Pauli 
matrices, and the remaining factors equal to the identity I2. The CCD is the decomposition 

u(2 JV ) = iX Q ffi il e (8) 

of 11(2^). The Lie subgroup e %1 ° associated to the subalgebra il is a subgroup of U(2 N ) containing 
all the local transformations. For each L £ U(2 N ) the decomposition ([2|) holds with K £ e lX ° and 
P = e p with P £ iT e . The factor K and in particular, any local transformation, does not modify 
the iV-qubit concurrence [1] . Such a decomposition is of type All if N is odd, and of type AI if N 
is even. 

The Odd-Even Decomposition (OED) was introduced in [4] as a generalization of the CCD to 
multipartite systems of arbitrary dimensions. The main idea is to construct a decomposition for 
the whole Lie algebra u(nin2 ■ ■ - tin) by combining decompositions for the Lie algebras associated 
to the single subsystems u(rij), j = l,...,N. This is based on the following observation for the 
CCD. When writing 

u(2) = span{icr x , ia y , ia z } © span{i I2}, 

we perform a (trivial) All decomposition of u(2), since su(2) = sp(l). In the CCD, we collect 
(modulo i) tensor products with an odd number of elements in the Lie algebra in X Q and tensor 
products with an even number of elements in I e . The OED [1] is obtained by applying this idea 
to general Lie algebras u(rtj), j = 1, 2, . . . , N. By writing 

u(rij) =K®V, 

with K. conjugate to so(rij) or sp(^f), and V — K , we obtain a decomposition of type AI or All, 
respectively. Denoting by a a generic element of ilC and by S a generic element of iV, we define 
X Q and X e to be the respective vector space spanned by tensor products of matrices of the type a 
and S with an odd or even number of a terms. The decomposition 

\x(nin 2 ■ ■ -n N ) — iX iX e (9) 

is a Cartan decomposition called the OED. The subspacc iX is the Lie subalgebra. This is a 
generalization of the CCD not only because it applies to systems of arbitrary dimensions, but 
also because, for every subsystem, we can perform different decompositions of type AI or AIL 
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The CCD is obtained as a special case of the OED © when all the subsystems are of dimension 
2 and a decomposition of type All is performed on each subsystem. Generalizing the result on 
the nature of the CCD decomposition, the OED decomposition is of type All if an odd number 
of All decompositions are performed. Otherwise, it is of type AI. As the CCD is related to the 
concurrence on N qubits, the OED has the same meaning for the generalized concurrences studied 
by Uhlmann in [lOj . 

We refer to [3] for a detailed discussion of the CCD and OED decompositions, and to [6] for 
the mathematical foundations of the Cartan decompositions. 

REMARK 2.1. The procedure described for the OED allows one great flexibility in the con- 
struction of various Cartan decompositions. Not only is one free to choose decompositions of type 
AI or All for each subsystem, but one can also choose among the different types of conjugate AI 
or All decompositions for each subsystem. This gives a method for the construction of an infinite 
number of decompositions in terms of tensor product matrices, even for the simplest case of N 
qubits. This flexibility is crucial in the construction of gradings for the Lie algebra u(nin2 • • -tin), 
and of recursive decompositions, as we shall see in the following two sections. 

We observe here that the procedure followed to construct the OED decomposition, applying 
decompositions of type AI and All, can be used with few changes to obtain an overall decompo- 
sition starting from decompositions of type AIII. More specifically, we perform decompositions of 
type AIII on each subsystem, and we collect in the respective subspaces il and iT e the linear 
combinations of tensor products with an odd or even number of factors in the subalgebra part 
(modulo i). We again consider the decomposition of u(n\n2 • ■ - tin) in @ but with il and il e 
defined in terms of type AIII decompositions. 

THEOREM 2.1. Consider the decomposition {P|) obtained with decompositions of type AIII as 
described above. This is a type AIII decomposition of the overall Lie algebra u{n\n2 ■ ■ -Tijv)- If 
N is odd, then il is the Lie subalgebra in the decomposition. If N is even, then iX e is the Lie 
subalgebra in the decomposition. 

Proof. The proof is by induction on N. If N = 1 the statement is obvious. Assume the statement 
is true for N — 1, and assume to be concrete that N is odd (exactly the same proof holds for N 
even). Denote by T J Q and I 3 e the respective spaces of matrices in iu{n\n2 • • -rij) that are linear 
combinations of an odd or even number of matrices in ifC, where K is the subalgebra of the 
AIII decomposition (possibly different for the different subsystems). Let us denote by /C the 
subalgebra of (block diagonal) matrices of the AIII decomposition on the last subsystem, and by 
V its orthogonal complement. We have 

i? = ( i?- 1 ® nc) e (i?- 1 ® iv) , i? = ( i?- 1 ®iv)®( i?- 1 ® tic 



By the inductive assumption, there exists a unitary matrix T in U{n\n2 ■ • ■ njsi-i) such that 
T^I^ I ~ 1 T is the same as the space of n\ni ■ ■ -njv-i x n\ri2 ■ ■ -n^-i Hermitian block-diagonal 
matrices, and T^I^~ l T is the same as the space of niri2 ■ ■ ■ nN-i x n\n2 ■ ■ ■ Jijv— l Hermitian 
block-antidiagonal matrices. Let T\ = (8 ln N , then the subspace T{T^Ti is spanned by all the 
matrices of the form 

A 0\ (C 0\ . / F\ ( G 
b)®{0 D)> and [-I* 0)®{-Gl 

The sizes of the matrices A, B, C, D, F, and G depend on the indices p, q of the two decompositions. 
Using the Corollary 4.3.10 of [13] one can construct a permutation similarity matrix T2 so that the 
subspace (TiT 2 ) t X ( ^ (TiT 2 ) is spanned by all the matrices of the form 

(C®A \ 

rpi ( A \ V - D ® A 

l2 \0 B)®{0 D) l2 ~ C®B 

\ D®B) 



G 



and by all the matrices of the form 



J 2 



F 
-Ft 



G 



F 2 = 



-Gt ' ' 

Finally, the conjugation P — > T3PT3, where T3 has the form 



/ 


-G®F+ 
\Gt ® Ft 



G®F\ 

-Gt ® F 



/ 



F 3 = 



/l 0\ 

1 

10 

\0 1 0/ 



with identity matrices 1 of appropriate dimensions, transforms the subspace (TiTsyig (F1F2) into 
the standard block diagonal form |(7J| of the type AIII decomposition. Therefore the subspace 



¥ into the standard 
□ 



is of form ([7]) where R := T1T2T3. It can be verified that R also transforms F; 
block anti-diagonal form of the AIII decomposition. 

REMARK 2.2. The indexes p and q of the resulting AIII decompositions of Theorem 12.11 are 
P = nin2 ■ ■ ■ ni-ipini+i ■ ■ ■ tin and q = n\ni ■ ■ ■ ni-iqini+i • ■ • tin, where I may refer to any of 
the subsystems I = 1,...,N, and pi and qi are the indices of the AIII decomposition of the l-th 
system. The theorem also indicates the inductive construction of the matrix conjugation which 
maps the AIII decomposition into the standard form. This is of interest for practical computation 
of the decomposition, as most of the existing numerical algorithms refer to the standard form (|6|), 
(0. Note that decompositions constructed by mixing AI or All decompositions with AIII- type 
decompositions do not give rise to Cartan decompositions. 



2.3 Recursive decompositions 

A recursive procedure to decompose unitary evolutions into local and entangling factors for the 
case of N qubits was introduced by N. Khaneja and S. Glaser in [7]. One starts with the Cartan 
decomposition 

su(2 JV )=/C®F (10) 

of type AIII, where 

K = span{l 2 <g> A , a z <g> B \ A G su(2 N - 1 ), B G u(2 Ar " 1 )} , 

V = span{cr x , y ®C\C<E u(2 JV_1 )}. l " ! 

This allows one to write each special unitary evolution L G SU(2 N ) as L = K1AK2, where 
K\, K2 £ e K and A G e A 7 with A the Cartan subalgebra contained in V . The subalgebra K. is the 
direct sum of spanji^ ® l 2 w-i} and two copies of su(2 Ar_1 ) which form a semisimple Lie algebra. 
Thus we can again apply Cartan's Theorem to further factorize K\,K2 G e . This is obtained 
through the decomposition K = K' ® V, with 

K! = span{l 2 «)A| Ae su(2 Ar ~ 1 )} , V' = span{cr z ® B \ B G u(2 jy - 1 )} , (12) 

to decompose each Fi ,F 2 G e K , thereby refining the decomposition of L. The key observation is 
that IC' and su(2 JV_1 ) are isomorphic, hence the procedure can be repeated by replacing N with 
N-l. 
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Another recursive procedure to decompose unitary evolutions was introduced by D. D'Alessandro 
and R. Romano in [5]. Such a decomposition applies to bipartite systems of arbitrary dimensions 
In the first step, one starts with an OED decomposition using AI types of decomposition on both 
subsystems, so that uiniriz) is decomposed as in Q, with 

1 := span{cr <g> S, S <g> a} , (13) 

conjugate to 50(774712). As 50(77477,2) is also semisimple, one then introduces a Cartan decomposition 
of 1 a by separating block diagonal and anti-diagonal elements (for two arbitrary indices) in the 
factors of the basis of T a . In particular, one writes 

do = K®V, 

where 

K := span{ia D (g> S D ,iS D <g> a D ,ia A <g> S A , iS A ® a A } 

and 

V := sp&n{ia D ® S A , iS D <g> icr A <g> S D , iS A <g> , 

the superscripts A and D standing for block-antidiagonal and block-diagonal respectively. The 
Lie algebra K is isomorphic to the semisimple direct sum so(r) so(f) with r + / = ni?i2. One 
decomposes /C as 

with /C' := span{icr D (g> S D 1 iS D <g> a D }, and V := spa,n{icr A ® S A , (8 cr A }. The Lie subalgebra 
/C' is isomorphic to the direct sum of four subalgebras So(n) 0so(r 2 ) 0so(r 3 ) ©so(r4). Each of the 
summands is spanned by tensor products of the type in (fT3")) with matrices cr and 5, where only 
one sub-block is different from zero. One then iterates the procedure. We refer to [5] for details. 



3 Lie algebra grading and recursive Cartan decompositions 

In this section, we give the definition of a grading of a Lie algebra, and relate a recursive Lie 
algebra decomposition to a grading. Our goal is to cast recursive decompositions of the unitary 
group into a general framework. In fact, in the following section we will show that known recursive 
decompositions, such as those of Khaneja and Glaser [7] and DAlessandro and Romano [S] reviewed 
in the previous section, can be obtained from an appropriate grading. We shall also see in the next 
section how new decompositions can be generated with the procedure described here. 

DEFINITION 3.1. Let C be a Lie algebra, and let M be an index set which has the structure 
of an additive semigroup. A direct sum decomposition 

C = (J) Ci 

is called an M-grading of £ if the subspaces Ci and Cj satisfy the commutation relation [d,Cj] C 
Ci+j for all i,j £ AI . 

In the special case where M is a monoid, that is, a semigroup with an identity element 0, the 
subspace Cq is a Lie subalgebra, since it satisfies the commutation relation [Cq, Cq] C Co . 

EXAMPLE 3.1. Consider the special linear Lie algebra sl(2) of 2x2 traceless matrices spanned 
by 

1\ _ fl \ _ (0 
o o) ' h ~ \0 -l) ' V ~ \ l 

with the commutation relations [h, x] — 2x , [x, y] = h, [h,y] = — 2y . Let M = {—1,0,1}. Then 
M becomes a monoid with addition given by the following table: 



1 Extensions to the general multipartite case can be obtained at the price of some notational complexity. 
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(M,+) 





- 1 1 








-1 1 


-1 


- 1 





1 


1 






The choice of £-± = spanja;}, Cq = span{/i} ; and C± = spanjy} makes sl(2) into an M-graded 
Lie algebra. 

A fundamental observation for what follows is that a Cartan decomposition ([1]) defines a Z2- 
grading of the Lie algebra £ with JC := Co and V ■= £i- 

As we have seen above, for a Lie algebra £, there are many Cartan decompositions. The 
following proposition shows that p Cartan decompositions give a Zj-grading for general p. 

PROPOSITION 3.1. Consider p Cartan decompositions £ = Cq © £{ for j = 1, . . . ,p. Define 

C klk2 ... kp := p| 4. (14) 

J=1,2,-.P 

for kj £ Z2 and j = 1, . . . ,p. TTien f/ie vector space decomposition 

£ = ®£ klk2 ... kp (15) 

forms a 1T 2 - grading of £. 

Proof. The proof is by induction on p. The claim is true for p = 1. Assume the claim is true for 
P — 1. Let ^4 £ £ kl k 2 ... kp and _B S £1^2. ..i p where € Z2 with 1 < i < p. Then it follows 

that A e £ klk2 ... kp _ 1 , A e £ kp , B e A 1 i 2 ..j p _ 1 , and B e £i p . Since £ is both Zf _1 -graded and 
Z2-graded, we have 

[A,B] e £(ki+h){k 2 +i 2 )-(k p -i+i p -i)> [A,B] € £( kp+ i p)l 

and therefore 

[A, B] e £(fci+; 1 )(fc2+*2)...(fc P -i+z f ,-i) n £( kp+ i p ) = C(ki+h)(k a +h)...(k p +ip), 
which implies that 

[^fcifca.-kp) Aii 2 ...i p ] C £(fe 1 +i 1 )(fe 2 +; 2 )...(fe t ,+/ p )- 
In conclusion, ([15]) is a Zf-grading of £. □ 

Thus a Cartan decomposition of a Lie algebra is a Z2-grading. A combination of p Cartan 
decompositions gives a Zj-grading. In order to cast a recursive decomposition in the framework of 
Lie algebra gradings, we give the following definition. 

DEFINITION 3.2. A recursive decomposition of a Lie algebra £ consists of two sequences of 
subspaces of £, 

So := {£o, £00, £000, ...,£op}, and Si := {£i,£oi,£ooi, £oj , - 1 i} 5 
both of length p, such that 

is a Cartan decomposition of £ j for each j — 0, ...,p— 1. That is, 

{£ j+i , £oi+ 1 ] Q £cp +i j [-£cp +i i Api] £ Apt) [A^ii Atfi] £ £o3+ 1 - 

-ffere we focroe se£ £o° : = £ an d £o°i := A- 
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Once one has a recursive decomposition of a Lie algebra £, in the sense of the above definition, 
one can obtain a decomposition of the connected Lie group e associated to C. This is obtained 
by repeated use of the Cartan decomposition theorem. Assume that C is semisimple, and that all 
of the C j , j = 1, ....,p— 1 are also semisimple. One first writes each element X of e c as 

X = K 1 AK 2 , 

where K\ and K% are in e Co , while A belongs to the connected Lie group corresponding to a 
maximal Abelian subalgebra contained in C\. Then one applies the Cartan decomposition of Co in 
order to decompose K\ and K 2 , and so on. The resulting decomposition contains several factors. 
A Zj-grading of C induces a recursive decomposition of C of length p: 

PROPOSITION 3.2. Consider a TF^-grading C = of C. Then the sequences So '■= 

{C k}, and Si := {C k-ii}, defined by 

A)* '•= @"R-a k ,j h+ i,...,j P and £o fc - 1 i : = R o'-'ij i+1 ,.,j I , 
for k = 1, ...,p, yield a recursive decomposition of C of length p. 



REMARK 3.1. Given a recursive decomposition sequence as in Definition [321 the semisimplicity 
of the subalgebras C k (for k = 0, ...,p) has to be verified independently. Even in the main 
case considered here, where the recursive decomposition sequence is obtained from a Lie algebra 
grading by means of combined Cartan decompositions as in Proposition 13. 11 semisimplicity is not 
guaranteed. For example, by combining type AI and All decompositions of su(4) in the standard 
basis, one obtains Cqq = sp(2) Hso(4). This is not semisimple, having an element which commutes 
with the whole Lie algebra. If a Lie algebra is the direct sum of a semisimple Lie algebra and 
an Abelian ideal, the Cartan decomposition theorem can be extended in the same fashion as we 
extended decompositions of su(n) to decompositions of u(n) in Section [27X1 



4 A scheme for recursive decompositions of U(n) 
4.1 Special cases 

We now show that the recursive decompositions of Khaneja and Glaser [J and D'Alessandro and 
Romano [5], summarized in Section [2.31 form a special case of the above procedure. In particular, 
they are induced by an appropriate grading. 

Let us start with the decomposition of Khaneja and Glaser [7]. We construct a Lie algebra 
grading of su(2 ) using the prescription of Proposition 13. II Consider p Cartan decompositions 

su(2 N ) = L\ © C{ (16) 

of su(2 N ) for j — 1, ...,p, all of type AIII, where C\ and C\ are respectively equal to and C\ in 
(jTTJ) . Now Cq and £| are defined in the same way as C\ and C\, except for the fact that a x and 
er z are interchanged. Such a decomposition is conjugate to the standard type AIII decomposition, 
the conjugation having the form A — > T ® 1 2 n-iAT' eg) 1 2 n-i with T the 2x2 matrix which 
diagonalizes cr^H The summands Cq and C\ are given by 

Ll = span{A ®l 2 ®C,B®a z ®D\A,B e u(2), C, D e n{2 N - 2 ) 1 Tr(A $ C) = 0} , 

C\ = span{£ ® a x . y ®F\Ee u(2) , F e u(2 Ar " 2 )} . ^ 17 ^ 



2 There is nothing special about a x here. One could have chosen u y instead. 
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This decomposition is again conjugate to the standard type AIII decomposition under the permu- 
tation which exchanges the first and second positions. The decomposition Cq © C\ is the same as 
£q © C\, except for the fact that the roles of a x and a z are exchanged. The summands C\ and 
C\ are defined analogously to C\ and C\, using the third position in place of the second. The 
same holds for Cq and C\ , defined as Cq and C\. In this fashion, one can define p = 2N — 1 
decomposition^! and therefore a Z 2 -grading of 511(2^). The corresponding pair of sequences giving 
the recursive decomposition according to Proposition 13 . 21 is 

Cq , Ci, same as in (p~T|) . 

C ao = span{l 2 ®A \A G su(2 N - 1 )} , 

An = span{cr z ®B\B € u(2 Ar " 1 )} , 



(18) 



£ooo 


= span{l 2 <g> 1 2 <X 


A, 1 2 (g)tr z >giC\Ae su(2 N - 2 ),C e u(2 N - 2 )} , 


£001 


= span{l 2 ©ctx,j/ 


Z>F\,Feu(2 N ~ 2 )}, 


Cq2N- 


-3 = span{l 2 -iv-i ( 


E) A , 1 2 n-2 ®a z ®c\Ae su(2), C G u(2)} , 


Cq2N- 


-<*! = span{l 2 «-2 


®a x , y ® F\,F €u(2)}, 


Cq2N- 


-2 = span{l 2 N-i ( 


E) A,A £ su(2)} , 


Cq2N- 


-si = span{l 2 w-2 


<g) a* <g> F | , F G u(2)} , 


Cq2N- 


-i = span{l 2 «-i ( 




Cq2N- 


-2i = span{l 2 «-i 





This sequence of subspaces is the one corresponding to the Khaneja-Glaser decomposition. Notice 
in particular that each Lie subalgebra £ k (for k = 1, 2N — 1) is either semisimple, or the sum 
of a semisimple Lie algebra and an Abelian (in fact one-dimensional) subalgebra of elements which 
commute with the whole Lie algebra. Thus the Cartan Decomposition Theorem applies in each 
case (cf. Remark |3"7T1) . 



In order to obtain the recursive decomposition corresponding to the decomposition of D'Alessandro 
and Romano [5], one constructs a grading by combining three types of decomposition: 

1) An OED decomposition with a type AI decomposition on each system; 

2) OED decompositions constructed using type AIII decompositions on each factor as in The- 
orem [23 

3) Type AIII decompositions in the standard form (separating block diagonal and block antidi- 
agonal matrices). 

In particular, let Cq — iT and C\ = iT^-, where X Q is defined in (fl"3"|) . The summands L\ and C\ 
are the respective subspaces il e and il a , referred to in Theorem 12. II The summands Cq and C\ 
are the JC and V subspaces of a type AIII decomposition in standard coordinates, with p and q 
given by p = p\P2 + Pi 92 and q = q\P2 + <Zi<Z2- Here {pi, qi} and {p 2 , g 2 } are the indices for the 
type AIII decompositions used for Cq and C\. The summands Cq, Cf, Cq and C\ are constructed 
analogously to Cq, C\, Cq, and C\ respectively, with different indices {pi,qi} and {p 2 ,g 2 }. The 
same holds for Cq, C\, Cq, and C\, and so on. Each time, the indices {p\,q\\ and {p 2 ,g 2 } are 
changed, differing from the previous ones in order to avoid repetition of decompositions. With 
these decompositions, one can define a grading, and therefore a recursive decomposition. This 
decomposition corresponds to the one in [5]. 



3 We stop at p = 2N — 1 because Cqp is {0}. 
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4.2 Construction of new recursive decompositions 

It follows from the previous discussion that many recursive decompositions of u(n) (or su(n)) 
and therefore of U(n) (or SU(n)), can be obtained. Once one has a certain number p of Cartan 
decompositions, then a Z^— grading and therefore a recursive decomposition can be obtained. We 
have seen that known recursive decompositions are a special case of this general procedure. Cartan 
decompositions can be obtained for example by taking one type of decomposition, e.g., AI, and then 
use various conjugations. When dealing with multipartite systems, it is convenient to have Cartan 
decompositions given in terms of tensor products of matrices as the CCD and OED described in 
subsection 12.21 

As an example we construct a new recursive decomposition of evolutions on N q-bits here. We 
consider the following 2N decompositions on u(2 N ). 



1) A CCD decomposition so that 

Cl=il?, C\=il», (19) 

where (I^)is the same as I (I e ) in ([5]) with the superscript N denoting the number of positions 
considered; 



2) An OED decomposition with all 'local' decompositions of type All except the one on the 
N— th term which is of type AI and of the form 

u(2) = span{ic z } ® span{il2, icr x , icr y }; 

For the resulting decomposition, we have 

C 2 a =span{i2 e Ar ~ 1 <g> a z , il^' 1 ® {<r x , a y , 1 2 }} , 
Cj = span{i2",f _1 <g> a z , iT^ -1 ® {cr x , a y , 1 2 }}. 



ri. 
L 'l- 



3) Same as in 2) but with o~ z and o~ x interchanged to define Cq and L\ ; 

4) Same as in 2) but with the N-th position replaced by the N — 1-th position to define L\ and 



5) Same as in 4) with a x and a z interchanged; 



6) — * 2N-1) ... and so on moving toward the first position, alternating decompositions as in 2) 
and decompositions as in 3); 

2N) Same as in 2) with the first position replacing the N-th one. 
EXAMPLE 4.1. In the case N = 3 we have, with a denoting any possible Pauli matrix, 

£q = spanjier ® 1 2 <g> I2, i I2 ®v ® 1 2 , i 1 2 ® I2 ia ® a <£> a}, 
C\ = span{icr ® a ® 1 2 , i I2 ®v ® a, ia <g> J 2 i 1 2 ® I2 ® I2}, 
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£o = span{«cr <g 1 2 ®{l 2 ,a x ,a y }, i l 2 ®a ® {i 2 , cr x , cr,,}, i<r ® a ® <j z ,i 1 2 ® 1 2 ®cr z }, 
£1 = span{icr <g> cr (g> {l 2 ,a x ,a y }, i 1 2 <8> l 2 ®{l 2 ,a x ,a y }, zcrig i 2 <gcr z , i 1 2 ®a ® <r z }, 

Cq = span{icr <g i 2 <g{-*2, cr z , cr y }, i <g {l 2 ,a z ,a y }, icr <%> a <g> a x ,i l 2 <g> 1 2 ®cr x }, 

C\ = span{icr <g cr ® {i 2 , Cz, i 1 2 ® 1 2 <g{i 2 , C2, c y }, icr <g> i 2 ®cr x , i l 2 ®a ® a x }, 

£q = span{i{l 2 , cr^, cr,,} ® it (g) 1 2 , «{i 2 , cr K , cr,J <g 1 2 <gcr, icr z <g cr <g cr, icr z (g 1 2 (g 1 2 }, 
£® = span{i{ 1 2 , <T X1 a y } (g cr <g cr, i{ i 2 , cr^, cr,,} (g J 2 (g 1 2 , ia z (g cr (g 1 2 , icr z (g J 2 (gcr}. 

In the general case, with the decompositions u(2 N ) := C 3 Q © C{, j = 1, ...,2N one constructs 
a grading as in Proposition 13.11 and a recursive decomposition according to Proposition 13.21 The 
sequences of subspaces associated to the latter are given, for k = 0, . . . , N — lp 

£ Q 2k+i = span{zl ( f r_fc <g l 2 fe}, 

C Q 2k 1 = spanlil^"* 1 <g a z <g l 2 k-i}, 

£ 2 fc+ 2 = span{iJ^" fc_1 (g l 2 *+i, il^^ 1 <g cr z (g l 2fc }, 

£o 2fc +!i = span{il ( f r_fe_1 (g {a x ,a y } (g l 2 fc}. 

In order to apply this recursive decomposition for the recursive decomposition of the Lie group 
U(2 N ) we make the following two remarks. 

REMARK 4.1. The Lie subalgebra £ 2fc + 1 = span{i2^~ fe (gl 2 iv-fc}, < k < N—l, is isomorphic 
to iT^~ k which is conjugate to so(2 N ~ k ) or sp(2 N ~ k ~ 1 ) according to whether N — k is even or 
odd, respectively. Thus, in every case the Lie algebra is semisimple. On the other hand, the Lie 
subalgebra £ 2fc = span-fiX^ <g l 2 k, iX^~ k <g <j z (g l 2 f=-i} is isomorphic to u(2 N ~ k ), and the 
isomorphism is given by the map 

A(g)l 2 ki — >A, B <g> a z ® l 2 fc-i ' — >B, (21) 

where A e iT^~ h and B G iTf - *. This is a the direct sum of a semisimple Lie algebra and a 
one dimensional subspace whose elements all commute with the elements of the Lie algebra. In all 
cases, the Cartan decomposition theorem applies. 

REMARK 4.2. In applying Cartan theorem to obtain a KAK decomposition as in ([3]) we need 
to identify the rank and a Cartan subalgebra at each step. The decomposition 

£o 2fe = £o 2fc + 1 ©£o 2fe l; 

with k — 0, N — lisa decomposition of type AI or All (modulo the isomorphism in (|2"Tj) ) of 
u(2 N ~ k ), according whether N — k is even or odd, respectively. In the AI case the rank is 2 N ~ k . 
A maximal Abelian subalgebra is spanned by the subspace 

'Haii — spanfiH w-t <g> a z ® l 2 fe-i|. 

2 

Here we have used the following notation 

H := span{cra; (g> a x , a y (g a y , a z (g a z , 1 2 <g) 1 2 }, (22) 

4 If the factors on the left occupy all the N positions in the tensor products, the factors on the right do not 
appear. 
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and Hi denotes the set obtained by tensor products of I elements of 7i, that is, Tii = 7i ® • • • ® 7i, 
I timesH In the odd, All, case, the rank is 2 N ~ k ~~ 1 . A Cartan subalgebra in this case is given by 

% aii '■= sp&nUH N-k-i (g) 1 2 ® u z ® l 2 fc-i } • 

2 

The decomposition 

with k — 0, TV— 1, is a decomposition of so(2 A '~ fe ) or sp(2 N ~ x ) according to whether N — k is 
even or odd. In the first case, it is a decomposition of type Dili (we refer to [5] for decompositions 
of Lie algebras different from u(n)) which has rank 2 N ~ k ~ 2 . The Cartan subalgebra can be taken 
equal to 

"Hdiii '•= span{ffl w-fc-2 ® i 2 ^ ® 1 2 *}. 

In the second case, it is a decomposition of type CI and the associated rank is 2 Ar_fe ~ 1 . The Cartan 
subalgebra can be taken equal to 

T~tci '■= span{ffl iY-t-i ® <r x (g> l 2 fc}- 



5 An Example of Computation 

In this section, we use an example to discuss some of the computational issues arising in recursive 
decompositions. In particular, we focus on the application of the recursive procedure described in 
the previous section to a generalized SWAP operator X sw <G U(8). In the tensor product basis, 
the action of X sw is defined by 

X sw : \i) <g> \j) ® \k) i ^ |i) ® |*) ® |t) , 

where i,j,k = 0,1, refers to an orthonormal basis { |0), |1)} of the Hubert space of each of three two 
level systems. X sw is the cyclic left shift operator acting on three qubits. The matrix representation 
of this operator is given by 



X, 



1 
00100000 



00001000 
00000010 
01000000 
00010000 
00000100 
\0 1/ 



Our goal is to factorize X sw in terms of elementary matrices. Using the Cartan decomposition 
of the previous section (cf. example 14. ip . one can construct a grading and therefore obtain a 
recursive decomposition of u(8). Modulo isomorphisms, the sequences characterizing the recursive 
decomposition are given by 

So = {sp(4),u(4), S o(4),u(2), S p(l),u(l)}, Si - { S p(4) ± ,u(4) ± , S0 (4) ± ,u(2) ± , S p(l) ± ,u(l) ± }. 

(23) 

Most of the algorithms for the computation of decompositions of the unitary group are given in 
standard coordinates. To transform the problem into standard coordinates, one uses an orthogonal 



5 Using the fact that H is a commuting set and induction on I along with the formula 
[K ® L, M ® N] = [K, M] ® {L. N} + {K, M} ® [L, N] , 



it is easy to see that Ti.[ is also a commuting set. 
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change of basis. According to [T] the associated matrix is given by 



F = 



1 

71 



1 










,1 






1 




1 



10 




1 








-1 





-1 






1 








-1 




1 










1 

-1 



I) 



(24) 



This matrix is referred to as the finagler. After this change of coordinates, X sw takes the form 
t y jp „,uv, v -u®XL, where 



X sw — F X SW F, with X sw 



a o o o\ 

x , I 1 

sw I 1 

\0 1 0/ 



To perform the decomposition, we follow the sequence of subspaces in ([23]) . The first step is to 
compute the decomposition of X sw induced by the Cartan pair (sp(4), sp(4)- L ) of u(8). It can be 
verified that X sw is symplectic, i.e., X sw £ Sp(A), therefore its decomposition is trivial. Moreover, 
X sw is contained in the image of 17(4) embedded into 5p(4)0 and represented by X' sw in J7(4). 
Indeed, X' sw is not only unitary but orthogonal, i.e., X' sw G 50(4). Hence the decompositions 
induced by the Cartan pairs (u(4), u(4)- L ), and (so(4),so(4)- L ) are also trivial. The computational 
problem is now to find the decomposition X' sw = K[A'K' 2 induced by the Cartan pair (u(2), u(2)- L ) 
so that K[ and K' 2 are contained the image of 17(2) embedding into 50(4), and A' is the exponential 
of an element of the suitable Cartan subalgebra, i.e., A 1 = (^ JLx ) . Let us partition X' sw into 
2x2 blocks, i.e., 

An A12 
A21 A22. 



A' 



Choose K 2 = I4. Then X' sw decomposes as 



A' 



.4 

-B 





E- 1 



(25) 



where A + iB G C/(2). This equation is equivalent to two matrix equations 
An - iX 2 i = {A + iB)E, A 22 + tX 12 = {A + iB)E~ 1 , 

which implies that 



E 2 = (A 22 + i Ai 2 )- 1 (An - iX 21 ) = (J ^ 



Once E is determined from the last equation, we obtain A and B using (f25|) so that 

(1 -1 ON 
110 
11 
0-11/ 





1 














1 










1 






1 









(26) 



In the next step, we decompose K[ using the Cartan pairs (sp(l), sp(l) ± ) and (u(l), u(l)- 1 -). 
As final result, we obtain 

X sw — LiL 2 L 3 L4, (27) 



6 An embedding of U(n) into Sp(n) or SO(2n) is given by the map U + iV <— > ( _y Yj ) where U and V are real 
matrices. 
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where 



Li = ;4=(l8 -i I2 ®v y <8> I2) , L 2 = 7^ +* I2 «) ct z ) , 
A 3 = -^(ls+i l 4 ®cr y ) , L 4 = 1 2 ®A', 

where A' is defined in ([25)) . We map X stu in (|27[) back to the tensor product basis to write 

X sw — L1L2L3L4, (28) 

where Lk = FLkF T , 1 < k < 4, where -F is the finagler defined in (f2~4")l . Finally, we write all the 
factors in (I28P as exponentials of matrices in the tensor product basis to obtain 



e 4 



1 e 4 ' 



e 4 



(29) 



For sake of comparison, we factorize X sw using the decomposition of Khaneja and Glaser [7]. 
We have shown that this factorization corresponds to the sequences 

So = {£0, £o 2 , -£o 4 i-£o 5 }i S\ = {£1, £oi ; £o 2 i, £o 3 i,£o 4 i} 

(cf. the elements of Sq and Si in (fl~8|) with AT = 3). Recall that £0 = span{l ®A, ® A> | A e 
su(4), _B e u(4)}. We find convenient to choose the Cartan subalgebra as the span of matrices of 
type a x ® D with D diagonal for the Cartan pair {Cq,C\) of u(8). Therefore the corresponding 
decomposition of X sw is given by 

X SW =K 1 AK 2 , 

where Kj = diag(Kji, Kj 2 ), with Kjk, 1 < j, k < 2, 4 x 4 unitary and A — where Dj 

diagonal with — £>| = I4. Following the procedure described in [3] (section 8.2.3), we obtain 
the matrices 



(I 



-1 

\0 

/o 



0\ 

1 



-1 0/ 



ATl 



/ 1 Q\ 
0-10 

-10 

\ 1/ 



#21 
















— i 

















— i 


\o 





1 


0/ 



AT- 



22 



A = 



D 2 = 



/l 0\ 


\0 1/ 

The next step is to factorize both K\ and K 2 using the Cartan pair (£ 2 , £ i )■ Notice that C\ 
is given by 1 <X>su(4). Choosing the Cartan subalgebra as the span of matrices of type a z ®D, with 
diagonal D, induces the decomposition 



-1 0\ 

i 

-i 

\0 1/ 



/o 0\ 

t 

-i 

\o oy 























A: 



-;2 





L j2 



^2 7 v - ~ JV v ~ ~j 

where Lji,Lj 2 £ SU(4) and A,- diagonal. In order to achieve this decomposition, we set 



(30) 



Kji 







a: 



j 2 



a: 
K 



p 
pt 



with unitary AT and P to obtain two matrix equations Kji = KP and Kj 2 — KP^ . Then it follows 
that P 2 = Kj 2 Kji. We diagonalize P 2 with a unitary matrix U to write P 2 = UXU\ and we 
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choose D = A' with det(D) — 1 so that P = UDW. Once P is determined, K can be found from 
the matrix equation Kji = KP^ . Finally we choose Lji = KU, Lj2 = IP and Aj = D to obtain 
the desired decomposition ([501) . 

Applying this procedure, we obtain 



Lll ~ 71 



(i 








1\ 





— i 


1 













-1 


V» 




1 


0/ 



/-l 1 



Ll2 ~71 







1 1 



-1 -1 
-1 1 

j 



(i 

1 

1 

\0 -i 



— i 




L 2 



1 

71 



/-l -1 0\ 

0-11 

11 

\-l 1 0/ 



Ai =A 2 



/i 0\ 

-i 

10 

\0 1/ 



Similarly, we repeat the first two steps with the respective Cartan pairs (£ 3 , £o 2 i) an d (£o 4 , £-o 3 i) 
to decompose Ljk, with 1 < k < 2. Finally, writing all the factors as exponentials, we obtain the 
factorization 

X sw = KiAK 2 , (31) 



with 



Ki 



^ I4 ®<y z 1 IgXTzlgXTz 1 g(7 x (g|l 1 4 H 1 4 (g)<T x 1 ®<T Z ®CT Z 



' I4 ® CTj; 



A= e 

'■e 



^ 2 _ g^f 1 I4 ®<r„ 1 ®<r z <g>lg-if 1 ®o z ®Oz p &Z- 1 4 p ^f 1 gajglp^ 1 1 ®c x ig)cr a 



1 4 gXTzg^jf 1 ®o'z®c Z g J f : 1 (glcr^gjlg^- 



'e 4 e 

■ 1 (g)Cr x ®(T z ,. I? I4 ig><7 



"e 2 



where we used 1 to denote 1.2- 



6 Conclusions 

Grading of a Lie algebra, Cartan decompositions and recursive decompositions of a Lie group are 
interrelated ideas. From a set of p Cartan decompositions, one can naturally obtain a Zj— grading 
of a Lie algebra and a recursive decomposition of the associated Lie group. Known procedures 
for the recursive decomposition of the unitary group of quantum evolutions are special cases of 
this general scheme. When dealing with multipartite quantum systems, it is convenient if the 
decompositions used in the procedure are given in terms of tensor products of basis elements of the 
Lie algebras associated to the single subsystems. This is the case for the CCD decomposition on 
n-qubits and the OED decomposition in its various forms. In this way, the factors of each element 
of the group are exponentials of tensor products, and one can identify local operations as well as 
multi-body interactions. 

We have given a new recursive decomposition applying the general procedure, along with an 
example of computation (section [5]) . For this example , formulas ((29]) and (f3Tj) . which is obtained 
applying the results of [7], give different decompositions. In general, different recursive decom- 
positions of u(n) will result in different factorizations of U(n). The framework presented here 
gives a virtually unbounded number of alternatives to decompose U (n) and parametrize quantum 
evolutions. 
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